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ABSTRACT 

Highly obscured aetive galaetic nuclei (AGN) are eommon in nearby galaxies, but are difficult to observe beyond the loeal Universe, 
where they are expected to signifieantly contribute to the black hole aecretion rate density. Furthermore, Compton-thick (CT) absorbers 
(Nh ~ 10^"^ em“^) suppress even the hard X-ray (2-10 keV) AGN nuelear emission, and therefore the column density distribution above 
10^"^ cm“^ is largely unknown. We present the identifieation and multi-wavelength properties of a heavily obscured (Nr ~ 10^^ cm“^), 
intrinsically luminous (L 2 -io> 10"^"^ erg s“^) AGN at z = 0.353 in the COSMOS field. Several independent indieators, sueh as the 
shape of the X-ray speetrum, the deeomposition of the spectral energy distribution and X-ray/[NeV] and X-my/6jLim luminosity ratios, 
agree on the fact that the nuelear emission must be suppressed by a ^ 10^^ cm"^ column density. The host galaxy properties show 
that this highly obscured AGN is hosted in a massive star-forming galaxy, showing a barred morphology, which is known to correlate 
with the presenee of CT absorbers. Finally, asymmetrie and blueshifted eomponents in several optieal high-ionization emission lines 
indieate the presence of a galactic outflow, possibly driven by the intense AGN activity (L^oilLEdd - 0-3 - 0.5). Such highly obscured, 
highly aecreting AGN are intrinsieally very rare at low redshift, whereas they are expected to be mueh more common at the peak of 
the star formation and BH accretion history, at z ~ 2 - 3. We demonstrate that a fully multi-wavelength approach can recover a sizable 
sample of such peculiar sources in large and deep surveys sueh as COSMOS. 
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1. Introduction 


Highly obscured, Compton thick (CT, column density Nr ~ lO^'^ 
cm“^) AGN are common in the local Universe, representing 
~ 20 - 30% of the local AGN population (Maiolino et al. 1998; 
Risaliti et al. 1999; Burlon et al. 2011). Several AGN-galaxy 
co-evolutionary models propose that their high-redshift counter¬ 
parts trace a crucial step in the AGN and galaxy build-up, cor¬ 
responding to the most efficient phase of black hoe (BH) growth 
and star formation activity (e.g. Hopkins et al. 2008), possibly, 
but not always, triggered by major mergers (Draper & Ballan- 
tyne 2012). A large population of CT AGN is also required to 
reproduce the shape of the cosmic X-ray background (Comas- 
tri et al. 1995; Gilli et al. 2007). More recently, the existence 
of a large population of heavily CT (HCT, Nh= 10^^“^^ cm“^) 
AGN has been proposed (Comastri et al. 2015) to reconcile the 
increased normalization of the local super massive black hole 
(SMBH) mass density (Graham & Scott 2013; Kormendy & Ho 
2013) with the normalization expected from growing black holes 
integrated over the cosmic time assuming an accretion efficiency 
of 6 ~ 0.1. Identifying these extremely obscured sources is, by 
definition, very difficult, even in the local Universe. An example 
of this is represented by Arp220: after years of debate (Genzel et 
al. 1998; Iwasawa et al. 2001; Clements et al. 2002; Nardini et 
al. 2010), ALMA observations have finally shown that this pro¬ 
totypical hyper-luminous infrared (IR) galaxy indeed hosts an 
AGN with an estimated column density of Nh=0.6 - 1.8 X 10^^ 
cm“^ (Wilson et al. 2014; Scoville et al. 2015). 


The identification of CT AGN beyond the local Universe is 
even more challenging. Recent results on preliminary NuSTAR 
data (Harrison et al. 2013) were only able to place an upper limit 
(< 33%) on the fraction of CT quasar (Liq-aqucV > erg s“^) 
between z = 0.5 - 1 (Alexander et al. 2013). Moreover, given 
the lack of large, complete samples of distant AGN observed 
at 20 - 30 keV, the relative fraction of HCT and CT sources 
remains poorly constrained and usually is not explored at all in 
deep survey studies (Tozzi et al. 2006; Ueda et al. 2014; but see 
also Brightman et al. 2014 and Buchner et al. 2015). 

In this article we present the discovery of a bona-fide HCT 
candidate at z = 0.353 in the COSMOS field (Scoville et al. 
2007), together with its multi-wavelength properties. Through¬ 
out this paper we assume i7o = 70 km s“^ Mpc“\ Oa = 0.73, 
and Qm = 0.27, and a Chabrier (2003) initial mass function. 


2. Selection 

The extreme properties of the source XMMUJ095910.4-1-020732 
in the XMM-Newton catalogue (Cappelluti et al. 2007^1 XID- 
392 hereafter, were serendipitously discovered by comparing the 
intrinsic (absorption-corrected) 2-10 keV luminosity (L 2 - 10 ) ob¬ 
tained from the automated spectral fit of all the sources in the 
XMM-Newton catalogue (Lanzuisi et al. 2015, LI5 hereafter), 
to the bolometric luminosity (Lboi) computed by fitting the spec- 

^ Chandra ID: 669, Elvis et al. (2009); optical RA=09:59:10.32, 
DEC=-h02:07:32.3, z = 0.353, Brusa et al. (2010); Civano et al. (2012). 
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Fig. 1: Left: The intrinsic L 2 - 10 , corrected for absorption is plotted against the bolometric luminosity derived from the SED fitting. 
Orange (cyan) circles represent type 2 (type 1) AGN, while green circles represent galaxies. The average error is shown in the upper 
left corner. Source XID-392 is shown with a red diamond and relative error bars. The relation derived for type 2 (type 1) AGN in 
L12 is shown with the red (blue) dashed line. Right: HST-ACS I band 20 x 20" image of the source. The X-ray contours from the 
Chandra full-band (0.5-7 keV) image are superimposed in green. 


tral energy distribution (SED) for the Herschel- detected sources 
in the same field (Delvecchio et al. 2014, 2015). Figure 1 shows 
the distribution of the two quantities for all the sample of XMM- 
Newton sources detected by Herschel. 

The Lboi of Fig. 1 (left) represents the total intrinsic lumi¬ 
nosity of the accretion disk, derived from the AGN model of 
the best-fitting SED solution. The SED-fitting procedure uses a 
grid of torus models that were computed by solving the radia¬ 
tive transfer equation for a smooth dusty structure irradiated by 
the accretion disc emission (Fritz et al. 2006; Feltre et al. 2012). 
This LboI does not include the X-ray emission above ^ 1 keV, 
which is negligible in the total budget, however. The two quan¬ 
tities plotted in Fig. 1 are related by the X-ray bolometric cor¬ 
rection ^Boi (Tbo 1 = ^BoiX T 2 - 10 ). The blue and red curves in Fig. 
1 represent the LboI- dependent ^boI relations derived in Lusso et 
al. (2012, LI2 hereafter) for type 1 and type 2 AGN in the COS¬ 
MOS field, respectively. Strikingly, source XID-392 is the only 
source in the sample (out of 394) that is more than two orders of 
magnitudes away from these relations, and -1.5 order of mag¬ 
nitudes away from the closest observed point for a given L 2-10 
or LboI- 

We verified that the source identification is correct and that 
the IR/optical and X-ray photometry is not contaminated by 
nearby sources. Figure 1 (right) shows the HST-ACS 20 x 20" 
cut-out of the source (Koekemoer et al. 2007). The host galaxy 
is an isolated barred spiral galaxy seen nearly face-on, and the 
closest source is at a distance of more than 7", which means that 
the IR/optical photometry is not contaminated. The X-ray con¬ 
tours, derived from the Chandra full-band (0.5-7 keV) imag^l, 
are superimposed in green. The off-set between the peak of the 
X-ray emission and the centre of the host galaxy (0.1") is well 
within the Chandra absolute astrometric accuracy (0.6" of ra¬ 
dius for the 90% uncertainty circle). Furthermore, XID-392 is 
the only X-ray emitting source within the 20 x 20" area. There¬ 
fore, the remaining possibility is that the absorption-corrected 

^ The HST and Chandra cutouts are publicly available at 
http://irsa.ipac.caltech.edu/data/COSMOS/index_cutouts.html 


L 2-10 is severely underestimated or that the LboI is severely over¬ 
estimated. 


3. Multi-wavelength properties 

3.1. Photometry and SED 

Figure 2a shows the rest frame broad-band SED fit, performed 
following Berta et al. (2013). The source is the third-brightest 
source at 24pm in the XMM-Newton catalogue and among the 
10% of the brightest sources in the PACS bands. The best-fit 
template includes a prominent AGN torus component (plotted 
in red) that dominates the source emission in all the IRAC (3.6, 
4.5, 5.8, and 8 pm) and MIPS 24 pm bands (the AGN fraction 
between 5 and 40 pm is 90%), therefore, the estimate of Lboi 
appears to be robust. The galaxy is massive (Log(MH=)=11.4 M©), 
with a derived star formation rate (SFR), after subtracting the 
AGN contribution in the IR, of SFR= 22.2 M© yr“^ The host 
galaxy lies at the massive end of the main sequence (MS) of star¬ 
forming galaxies (e.g. Whitaker et al. 2012) at z=0.353, having 
a specific SFR (sSFR=SFR/M,) of sSFR= 9.3 x lO'^ Gyr'^. 

The source is also detected as a compact source at 1.4 GHz 
in the VLA COSMOS survey (Schinnerer et al. 2004), with a 
flux of Fiaghz - 0.424 + 0.025 mJy. Given the definition of 
the q parameter from Helou et al. (1985, q = Log(FIRI3.15 x 
lO^^Hz)/Sy(1.4GHz)), and a far-infrared (FIR, i.e. rest frame 
42.5-122.5 pm) flux of FIR= 1.1 x 10“^^ erg cm“^ s“^ of the 
star-forming component, XID-392 has ^ = 1.8 + 0.1. Therefore, 
it has a radio flux that is higher than the average q observed in lo¬ 
cal and high-redshift star-forming galaxies {qmean = 2.21 +0.18 
from Del Moro et al. 2013; qmean = 2.26 ± 0.08 from Magnelli 
et al. 2015). A significant fraction (~ 45%) of these radio-excess 
galaxies are thought to host low-to-moderate luminosity, or lu¬ 
minous dust-obscured AGN (Del Moro et al. 2013). 
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Fig. 2: a: Rest frame, broad-band SED of XID-392. In blue we show the stellar emission, in red the AGN torus emission, and 
in magenta the SF dust emission, b: Unfolded XMM-Newton (black) and Chandra (red) spectra of XID-392. The best-fit model 
(magenta) is composed of a thermal component (green) and an obscured torus template (blue), c: SDSS optical spectrum of XID- 
392. The most prominent emission lines are labelled. The grey area masks a region of bad-sky subtraction, d: [NeV]3346,3426 and 
[FeVII]6086 lines decomposition. The green curve represents the narrow component, while the blue curve represents the broadened 
asymmetric component, associated with the outflow. The red curve is the sum of the two components. The dotted lines represent the 
systemic redshift, estimated from the continuum and the stellar absorption lines. 


3.2. X-ray data 

The Chandra and XMM-Newton spectra of XID-392 were ex¬ 
tracted and fitted according to Lanzuisi et al. (2013) and Mainieri 
et al. (2007). Figure 2b shows the Chandra and XMM-Newton 
spectra with the best-fit model described below. Both spectra 
look peculiar, despite the limited photon statistics (64 net counts 
in XMM-Newton and 44 in Chandra). They show soft emission 
extending up to 2 keV that can be modelled with either a soft 
power law with photon index F ~ 3.7 or with thermal emission. 
The source is not significantly detected in the 2-10 keV band 
in the XMM-Newton datfl while the Chandra data, thanks to 
the deeper exposure and lower background, enables the detec¬ 
tion also in the hard band (with ~ 10 net counts in the 2-10 keV 
band). Despite the limited spectral quality, the shape of the hard 
part of the spectrum is clearly flatter than the typical AGN power 
law with F = 1.9 (Piconcelli et al. 2005): fitting the spectra with 
a simple power law in the 2-10 keV band alone yields a best-fit 
photon index of F = 0.68 + 0.45. 

Assuming that the hard component is produced by obscura¬ 
tion of the primary powerlaw, it can be modelled with the torus 
template from Brightman & Nandra (2011, consistent results are 
obtained with my torus of Murphy & Yaqoob 2009). This model 
gives a best-fit Nr value close to 10^^ cm“^. Given the limited 


^ The source was not selected as CT in LI5 because of the lack of hard 
detection in the XMM-Newton catalogue. 


photon statistics available in the 2-10 keV band, however, this 
value is poorly constrained. To estimate at least a lower limit on 
Nh, we fixed the normalization of the torus template to the nor¬ 
malization required for an intrinsic Log(L2-io)=44.25 erg s“\ 
that is, the intrinsic L 2-10 expected from LboI, using the ^boI of 
LI2. With this constraint, the resulting lower limit is Nr ~ 5 x 
10^"^ cm“^. The 2-10 keV observed (not corrected for absorption) 
luminosity is instead well constrained .1 +0.6 erg 

S"'fl 

All nearby CT AGN show a strong {EW ^ 1 keV) emission 
line at the rest frame energy of the Fe Ka line at 6.4 keV. This 
feature is absent from the spectrum of XID-392. However, the 
very limited number of counts available between 6 and 7 keV 
rest frame (about four net counts in total) only allowed us to 
estimate a loose upper limit for the equivalent width of the Fe 
Ka line of EW < 1.4(2.4) keV at 90% (3cr) confidence level. 
The upper limit on the intensity of the Fe Ka emission line is 
therefore fully consistent with the possibility XID-392 is a CT 
AGN. 

Given the SFR derived in Sec. 3.1, the expected 
0.5-2 and 2-10 keV band luminosities from SFR are 


^ The L 2-10 plotted in Fig. 1 was computed from the automated fit de¬ 
scribed in L15 with a single power law and F = 1.9 fixed. The observed 
2-10 keV luminosity was therefore overestimated and, when consider¬ 
ing the spectral model adopted in this work, the source is even more 
extreme in terms of distance from the L 2 - 10 -TB 01 relation. 
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Lo. 5 - 2 (SFR)~L 2 -io(SFR)~ 10"^^ erg s“^ according to the relation 
reported by Ranalli et al. (2003) or Lq.s-sCSFR)- 9 x 10"^^ erg 
s“^from the Mineo et al. (2014) relation. These values are about 
one order of magnitude lower than what was observed from the 
Chandra and XMM-Newton spectra of XID-392 (L(Obs)=1.3, 
0.5, l.bxlO"^^ erg s“Hn the 0.5-2, 2-10 and 0.5-8 keV band, re¬ 
spectively). Therefore the presence of a second source of X-ray 
photons (i.e. the obscured AGN) is required to explain the ob¬ 
served luminosity. 

4. Obscuration diagnostics 

What we described above shows that the X-ray data alone, al¬ 
though consistent with the presence of an AGN obscured by 
an HCT absorber, do not provide conclusive evidence. The rich 
multi-wavelength data set, both spectroscopic and photometric, 
that is available in the COSMOS field, allows us to explore the 
properties of this unique source in more detail. 

4.1. NeV diagnostics 

Figure 2c shows the SDSS DRIO (Ahn et al. 2013) spectrum 
of XID-392. Several strong emission lines are visible in both 
the blue and red part of the SDSS spectrum. Unfortunately, the 
lOIII] emission line, a commonly used indicator of AGN activ¬ 
ity, is masked in the SDSS spectrum because of problems in the 
sky subtraction. Although it is typically weaker than 10III]5007 
and suffers stronger dust extinction, the lNeV]3426 line is con¬ 
sidered an unambiguous sign of nuclear activity (Schmidt et al. 
1998; Mignoli et al. 2013) because high-energy photons (> 0.1 
keV) are required to produce this line. 

Gilli et al. (2010) introduced a diagnostic based on the ra¬ 
tio between the observed rest-frame 2-10 keV band and the 
lNeV]3426 luminosities that was calibrated on a sample of ob¬ 
scured and unobscured local Seyferts. Figure 3 (left) shows the 
L^^_\QlL(NeV) ratio for several samples of obscured and unob¬ 
scured sources (Vignali et al. 2014; Gilli et al. 2010; Young et 
al. 2009). The solid line shows the expected trend obtained by 
starting from the mean ratio {X/[NeV]} observed in unobscured 
objects and progressively obscuring the X-ray emission with in¬ 
creasing Nh (up to log(NH) = 25.5) while keeping the INeV] 
luminosity fixed, because it is produced in the more extended 
narrow line region. The cyan (grey) shaded regions correspond 
to ±lcr (+90%) around the {X/[NeV]) ratio (Gilli et al. 2010). 
Source XID-392 is shown as a red diamond and has indeed a 
very strong INeV] (EW = 43 A, Fa^^v = 1.33 x 10“^^ erg cm“^ 
s“^) and an extreme value of X/lNeV]= 1.44 ± 1.20. The Nr 
lower limit is derived from the X-ray spectral fit (Sect. 3.2). Re¬ 
markably, XID-392 has the second-lowest X/lNeV] ratio of the 
sample, and its value is one order of magnitude lower than the 
threshold defined in Gilli et al. (2010) to select CT sources. This 
confirms the exceptional nature of this source. 

4.2. Mid-IR diagnostics 

The AGN intrinsic L 2-10 (T^O and the Ljr re-emitted by the ob¬ 
scuring torus are known to follow a tight correlation over sev¬ 
eral orders of magnitudes (Lutz et al. 2004; Gandhi et al. 2009). 
Given that the observed L 2-10 is affected by obscuration, 

while the Ljr is largely independent of it, the selection of AGN 
with very low to Ljr ratios has been extensively used to 
identify CT sources. The vs Ls.gyum relation is linear for 
low-redshift Seyfert galaxies, up to L 5 . 8 yum= 10"^^ erg s“^(Lutz 


et al. 2004), while several works on luminous distant QSOs have 
shown that the relation tends to flatten above L 5 . 8 yum= 10'^'^ erg 
s“^ (Maiolino et al. 2007; Fiore et al. 2009; Lanzuisi et al. 2009, 
L09 hereafter). 

Figure 3 (right) shows the distribution of vs. L 5 . 8 yum for 
several samples of CT candidates from the CDFN (Alexander 
et al. 2008), CDFS (Georgantopoulos et al. 2013), X-SWIRE 
(L09), and from recent NuSTAR observations (Lansbury et al. 
2014; Stern et al. 2014). The range of redshifts covered by 
the different samples is wide, and even if source XID-392 is 
at the lowest end of the redshift distribution, it is intrinsically 
bright, with L 5 . 8 yum~ 10^^ erg s“^ This means that it is just at 
the intersection of the local low-luminosity and the high-z high- 
luminosity regime. We stress that the source is the only one be¬ 
low the Log(NH)=25 cm“^ line (dotted line in Fig. 3 right) if the 
relation described by Fiore et al. (2009) is considered, while it 
would be even more extreme in terms of vs. L 5 . 8 yum for the 
relation reported by Lutz et al. (2004). 

5. Discussion 

The different diagnostics discussed above strongly indicate that 
the nuclear emission of XID-392 must be obscured by an HCT 
absorber (Nr ~ 10^^ cm"^). The optical spectrum allow us to 
explore other physical properties of this unique source. An esti¬ 
mate of the SMBH mass (M^//) can be obtained from the broad 
component of the Ha line in the optical spectrum and the 5100A 
continuum luminosity or the L 2 - 10 , using the relations reported 
in Bongiorno et al. (2014, Eqs. 2 and 4). The Htr-INII] com¬ 
plex has a total flux of 5.5 x 10“^^ erg cm"^ s“\ while the broad 
component of the Ha line has a flux of 2.7 x 10“^^ erg cm“^ s“^ 
and a FWHM of 2700 km/s, with large uncertainties due to the 
multiple line decompostion. The 5100A continuum luminosity 
is 1.7 X 10"^^ erg s“^ The resulting Mrh are LogiMBu) = 8.18 
Mq and Log{MBH) = 8.35 M© from the two relations, respec¬ 
tively. Consistent results are obtained if we estimate Mbh from 
the host derived from the SED fitting, rescaled to the bulge 
mass using the mean bulge-to-total mass ratio of 0.2 reported by 
Kormendy & Ho (2013) for local barred galaxies. This value is 
further rescaled to the using the BH-to-bulge mass ratio of 
0.0023 reported in Marconi & Hunt (20030 The derived value 
is Log{MBH) ~ 8.1 Mq. 

Given the LboI reported in Sect. 2, the Eddington ratio is 
7 ^Boi/^£JJ= 0-3-0.5, that is, at the upper end of the distribution for 
X-ray selected AGN in deep surveys. The sSER, close to the MS 
value, and the face-on barred spiral morphology indicate, how¬ 
ever, that XID-392 has not been caught in a major or gas-rich 
merger event, but is instead an HCT in which the obscuration 
is likely to take place in a small-scale torus and not in the host 
or starburst regions. Indeed, the presence of a bar in the host is 
known to correlate with CT obscuration (Maiolino, Risaliti & 
Salvati 1999) and with AGN activity in general (Galloway et al. 
2015). 

The optical spectrum of XID-392 is complex, showing 
strong Mgll, INeV], tOII], INelll], INII], and Ha emission 
lines, as well as a strong 10III]3133A Bowen fluorescence line 
(Schachter et al. 1990) and the high-ionization lEeVII]6086A 
emission line. Eurthermore, almost all these emission lines show 
an asymmetric profile, and remarkably, both the INeV] and lEe- 
VII], having a similarly high ionization potential of ~ 125 eV. 

^ The BH-to-bulge mass ratio computed in Kormendy & Ho (2013) 
was derived for classical bulges and cannot be applied to XID-392, 
which is a barred spiral. 
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Fig. 3: Left: Rest frame 2-10 keV to [NeV] luminosity ratio as a function of Nr. Source XID-392 is shown as a red diamond. V14: 
Vignali et al. 2014, GIO: Gilli et al. 2010, Y09: Young et al. 2009. Right: Log(L^^0 vs LogCLs g/^m) for several CT candidate (L14: 
Lansbury et al. 2014, G13: Georgantopoulos et al. 2013, L09: Lanzuisi et al. 2009, Stern et al. 2014, AOS: Alexander et al. 2008, 
P15: Piconcelli et al. 2015). The green shaded area is the relation reported by Lutz et al. (2004) for a sample of low-redshift AGN. 
The black solid line is the relation for the high-redshift AGN described by Fiore et al. (2009), the long dashed (dotted) line the 
expected relation for a 10^^ cm“^ (10^^ cm“^) absorber, starting from the relation given by Fiore et al. 


Figure 2d shows a zoom in the [NeV] and [FeVII] region. The 
emission lines are fitted with a narrow component (green curve) 
plus a broadened component (blue curve). Both components are 
blueshifted with respect to the systemic velocity, estimated from 
the continuum and the stellar absorption lines, with a typical ve¬ 
locity offset of ~ 300 and ~ 800 km/s, respectively. The FWHM 
of the broadened component is in the range 1200-1800 km/s. 
These features are commonly associated with outflowing gas on 
kpc scales (Harrison et al. 2012, 2014), and similar features have 
been found in a handful of sources in the XMM-COSMOS cat¬ 
alogue, with similar accretion properties, that is, high L^oilLEdd 
(Brusa et al. 2010; Perna et al. 2015). The detailed characteri¬ 
zation of the properties of the putative outflow are the subject 
of ongoing observational effort, and the results will be presented 
in a forthcoming paper. We stress, however, that the exceptional 
amount of obscuration and the strong outflow signatures are pos¬ 
sibly related: if the AGN is in the short-lived phase in which 
powerful accretion (Lboi/^£^/j= 0.3-0.5) occurs in a surrounding 
dense gas cocoon (Hopkins et al. 2005, 2008), strong outflow 
episodes, driven by radiation pressure from the AGN, are ex¬ 
pected (Menci et al. 2008). 

The presence of extreme obscuration and strong outflow sig¬ 
natures in the same source does not require any particular geom¬ 
etry, except for the fact that the line of sight must intercept the 
obscuring material. A local example of such a configuration is 
NGC 1068, which is an HCT (Nh> 10^^ cm'^ Matt et al. 2004) 
and shows a fast (Ymax = 1400 km/s) outflow with an outflow 
half-opening angle 6out = 27° and an inclination angle of i = 9° 
almost perpendicular to the line of sight (Muller-Sanchez et al. 
2011). In a similar scenario that does not imply a wide opening- 
angle for the outflow, the blueshift observed in the emission lines 
is only the radial component of the real outflow velocity, and in¬ 
deed the velocities inferred from the blueshift should be consid¬ 
ered as lower limits. 

Detecting one such highly obscured, highly accreting source 
will not change our understanding on the demographics of HCT, 


of course. However, it is important to recognize that these 
sources do exist beyond the local Universe, and they are indeed 
expected to be much more common at the peak of star forma¬ 
tion or BH accretion history, at z ~ 2-3. More importantly, 
we have demonstrated that they are detectable, if a large multi¬ 
wavelength survey is combined with medium-deep X-ray obser¬ 
vations and we properly combine all the information simultane¬ 
ously. A systematic approach in this direction should be able to 
retrieve a non-negligible population of such HCT sources in a 
wide range of redshift (thanks to the positive K-correction that 
applies to obscured AGN in X-rays) if applied to a large area, 
multi-wavelength survey such as COSMOS: the X-ray back¬ 
ground synthesis model described by Gilli et al. (2007) predicts 
~ 2.2 deg“^ of these HCT sources (LogNH=25-26 cm“^), as¬ 
suming a flux limit of F 2 -io=l x 10“^^ erg cm“^ s“\ that is, 
the flux observed for a source like XID-392. Therefore the full 
COSMOS-Legacy survey (2 deg^, Civano et al. 2015a submit¬ 
ted) will deliver a small but valuable sample of such sources. 
We stress that a source like XID-392 is far below the detectabil¬ 
ity above 10 keV with Swift-EKL or INTEGRAL-IBIS, while 
NuSTAR would require a deep > 500 ks exposure, which is not 
feasible on a large-area survey like COSMOS. The source is not 
detected in the current ~ 100 ks exposure available in the field 
(Civano et al. 2015b submitted), which already required a to¬ 
tal of 3.2 Ms of NuSTAR observing time, therefore the multi¬ 
wavelength approach is the only feasible approach for this class 
of sources with the current instrumentation. 
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